iMapSplice: Alleviating Reference Bias Through Personalized RNA-seq Alignment

Genomic variants in both coding and non-coding sequences can have functionally important and sometimes deleterious effects on exon splicing of gene transcripts. For transcriptome profiling using RNA-seq, the accurate alignment of reads across exon junctions is a critical step. Existing algorithms that utilize a standard reference genome as a template sometimes have difficulty in mapping reads that carry genomic variants. These problems can lead to allelic ratio biases and the failure to detect splice variants created by splice site polymorphisms. To improve RNA-seq read alignment, we have developed a novel approach called iMapSplice that enables personalized mRNA transcriptome profiling. The algorithm makes use of personal genomic information and performs an unbiased alignment towards genome indices carrying both reference and alternative bases. Importantly, this breaks the dependency on reference genome splice site dinucleotide motifs and enables iMapSplice to discover personal splice junctions created through splice site polymorphisms. We report comparative analyses using a number of simulated and real datasets. Besides general improvements in read alignment and splice junction discovery, iMapSplice greatly alleviates allelic ratio biases and unravels many previously uncharacterized splice junctions created by splice site polymorphisms, with minimal overhead in computation time and storage. Software download URL: https://github.com/LiuBioinfo/iMapSplice

Cellular Proliferation of Equine Bone Marrow- and Adipose Tissue-Derived Mesenchymal Stem Cells Decline With Increasing Donor Age

Background: Bone marrow (BM)- and adipose tissue (AT)-derived mesenchymal stem cells (MSCs) are used increasingly for autologous cell therapy in equine practice to treat musculoskeletal and other injuries. Current recommendations often call for 10–100 million MSCs per treatment, necessitating the expansion of primary cells in culture prior to therapeutic use. Of concern, human and rodent studies have shown a decline of both MSC recovery from sampled tissue and in vitro proliferative capacity with increasing donor age. This may be problematic for applications of autologous cell-based therapies in the important equine demographic of older patients. Objectives: To investigate the effect of donor age on the cellular proliferation of equine BM- and AT-MSCs. Study Design: In vitro study. Methods: BM- and AT-MSCs and dermal fibroblasts (biological control) were harvested from horses in five different age groups (n = 4, N = 60); newborn (0 days), yearling (15–17 months), adult (5–8 years), middle-aged (12–18 years), and geriatric (≥ 22 years). Proliferation of the cells was tested using an EdU incorporation assay and steady state mRNA levels measured for targeted proliferation, aging, and senescence biomarkers. Results: The cellular proliferation of equine BM- and AT-MSCs declined significantly in the geriatric cohort relative to the younger age groups. Proliferation levels in the two MSC types were equally affected by donor age. Analysis of steady state mRNA levels showed an up-regulation in tumor suppressors, apoptotic genes, and multiple growth factors in MSCs from old horses, and a down-regulation of some pro-cycling genes with a few differences between cell types. Main Limitations: Potential age-dependent differences in cell function parameters relevant to cell-therapy application were not investigated. Conclusions: The cellular proliferation of equine BM- and AT-MSCs declined at advanced donor ages. High levels of in vitro proliferation were observed in both MSC types from horses in the age groups below 18 years of age

Discerning Novel Splice Junctions Derived from RNA-Seq Alignment: A Deep Learning Approach

Background: Exon splicing is a regulated cellular process in the transcription of protein-coding genes. Technological advancements and cost reductions in RNA sequencing have made quantitative and qualitative assessments of the transcriptome both possible and widely available. RNA-seq provides unprecedented resolution to identify gene structures and resolve the diversity of splicing variants. However, currently available ab initio aligners are vulnerable to spurious alignments due to random sequence matches and sample-reference genome discordance. As a consequence, a significant set of false positive exon junction predictions would be introduced, which will further confuse downstream analyses of splice variant discovery and abundance estimation. Results: In this work, we present a deep learning based splice junction sequence classifier, named DeepSplice, which employs convolutional neural networks to classify candidate splice junctions. We show (I) DeepSplice outperforms state-of-the-art methods for splice site classification when applied to the popular benchmark dataset HS3D, (II) DeepSplice shows high accuracy for splice junction classification with GENCODE annotation, and (III) the application of DeepSplice to classify putative splice junctions generated by Rail-RNA alignment of 21,504 human RNA-seq data significantly reduces 43 million candidates into around 3 million highly confident novel splice junctions. Conclusions: A model inferred from the sequences of annotated exon junctions that can then classify splice junctions derived from primary RNA-seq data has been implemented. The performance of the model was evaluated and compared through comprehensive benchmarking and testing, indicating a reliable performance and gross usability for classifying novel splice junctions derived from RNA-seq alignment

Analysis of equine protein-coding gene structure and expression by RNA-sequencing

RNA-sequencing (RNA-seq) data from eight equine tissue samples (34-day whole embryo, full term placental villous, adult testes, adult cerebellum, adult articular cartilage, adult LPS-stimulated articular cartilage, adult synovial membrane, and adult LPS-stimulated synovial membrane) were used to refine the structural annotation of protein-coding genes in the horse and for a preliminary assessment of tissue-specific expression patterns

Differential gene expression associated with postnatal equine articular cartilage maturation

<p>Abstract</p> <p>Background</p> <p>Articular cartilage undergoes an important maturation process from neonate to adult that is reflected by alterations in matrix protein organization and increased heterogeneity of chondrocyte morphology. In the horse, these changes are influenced by exercise during the first five months of postnatal life. Transcriptional profiling was used to evaluate changes in articular chondrocyte gene expression during postnatal growth and development.</p> <p>Methods</p> <p>Total RNA was isolated from the articular cartilage of neonatal (0–10 days) and adult (4–5 years) horses, subjected to one round of linear RNA amplification, and then applied to a 9,367-element equine-specific cDNA microarray. Comparisons were made with a dye-swap experimental design. Microarray results for selected genes (COL2A1, COMP, P4HA1, TGFB1, TGFBR3, TNC) were validated by quantitative polymerase chain reaction (qPCR).</p> <p>Results</p> <p>Fifty-six probe sets, which represent 45 gene products, were up-regulated (p < 0.01) in chondrocytes of neonatal articular cartilage relative to chondrocytes of adult articular cartilage. Conversely, 586 probe sets, which represent 499 gene products, were up-regulated (p < 0.01) in chondrocytes of adult articular cartilage relative to chondrocytes of neonatal articular cartilage. Collagens, matrix-modifying enzymes, and provisional matrix non-collagenous proteins were expressed at higher levels in the articular cartilage of newborn foals. Those genes with increased mRNA abundance in adult chondrocytes included leucine-rich small proteoglycans, matrix assembly, and cartilage maintenance proteins.</p> <p>Conclusion</p> <p>Differential expression of genes encoding matrix proteins and matrix-modifying enzymes between neonates and adults reflect a cellular maturation process in articular chondrocytes. Up-regulated transcripts in neonatal cartilage are consistent with growth and expansion of the articular surface. Expression patterns in mature articular cartilage indicate a transition from growth to homeostasis, and tissue function related to withstanding shear and weight-bearing stresses.</p

Functionally Annotating Regulatory Elements in the Equine Genome Using Histone Mark ChIP-Seq.

One of the primary aims of the Functional Annotation of ANimal Genomes (FAANG) initiative is to characterize tissue-specific regulation within animal genomes. To this end, we used chromatin immunoprecipitation followed by sequencing (ChIP-Seq) to map four histone modifications (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in eight prioritized tissues collected as part of the FAANG equine biobank from two thoroughbred mares. Data were generated according to optimized experimental parameters developed during quality control testing. To ensure that we obtained sufficient ChIP and successful peak-calling, data and peak-calls were assessed using six quality metrics, replicate comparisons, and site-specific evaluations. Tissue specificity was explored by identifying binding motifs within unique active regions, and motifs were further characterized by gene ontology (GO) and protein-protein interaction analyses. The histone marks identified in this study represent some of the first resources for tissue-specific regulation within the equine genome. As such, these publicly available annotation data can be used to advance equine studies investigating health, performance, reproduction, and other traits of economic interest in the horse

Analysis of unannotated equine transcripts identified by mRNA sequencing

Sequencing of equine mRNA (RNA-seq) identified 428 putative transcripts which do not map to any previously annotated or predicted horse genes. Most of these encode the equine homologs of known protein-coding genes described in other species, yet the potential exists to identify novel and perhaps equine-specific gene structures. A set of 36 transcripts were prioritized for further study by filtering for levels of expression (depth of RNA-seq read coverage), distance from annotated features in the equine genome, the number of putative exons, and patterns of gene expression between tissues. From these, four were selected for further investigation based on predicted open reading frames of greater than or equal to 50 amino acids and lack of detectable homology to known genes across species. Sanger sequencing of RT-PCR amplicons from additional equine samples confirmed expression and structural annotation of each transcript. Functional predictions were made by conserved domain searches. A single transcript, expressed in the cerebellum, contains a putative kruppel-associated box (KRAB) domain, suggesting a potential function associated with zinc finger proteins and transcriptional regulation. Overall levels of conserved synteny and sequence conservation across a 1MB region surrounding each transcript were approximately 73% compared to the human, canine, and bovine genomes; however, the four loci display some areas of low conservation and sequence inversion in regions that immediately flank these previously unannotated equine transcripts. Taken together, the evidence suggests that these four transcripts are likely to be equine-specific

Tissue Restricted Splice Junctions Originate Not Only from Tissue-Specific Gene Loci, but Gene Loci with a Broad Pattern of Expression

Cellular mechanisms that achieve protein diversity in eukaryotes are multifaceted, including transcriptional components such as RNA splicing. Through alternative splicing, a single protein-coding gene can generate multiple mRNA transcripts and protein isoforms, some of which are tissue-specific. We have conducted qualitative and quantitative analyses of the Bodymap 2.0 messenger RNA-sequencing data from 16 human tissue samples and identified 209,363 splice junctions. Of these, 22,231 (10.6%) were not previously annotated and 21,650 (10.3%) were expressed in a tissue-restricted pattern. Tissue-restricted alternative splicing was found to be widespread, with approximately 65% of expressed multi-exon genes containing at least one tissue-specific splice junction. Interestingly, we observed many tissue-specific splice junctions not only in genes expressed in one or a few tissues, but also from gene loci with a broad pattern of expression